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We discuss the occupation number correlations in an ultracold system of interacting fermionic 
atoms. For a system with a special energy-level distribution, viz. two multiply-degenerate levels, 
explicit expressions for the correlation functions are derived in a canonical approach using the exact 
ground state wavefunction of the reduced BCS Hamiltonian. We evaluate the correlators numerically 
for different interaction strength and find analytical expressions in some limiting cases. Due to the 
underlying fermionic nature of the pairs the occupations are predominantly anti-correlated and their 
statistics is a multinomial distribution. 



I. INTRODUCTION 



Ultracold fermionic gases have attracted considerable 
attention in theoretical and experimental physics re- 
cently. This has been intensified after the experimental 
successes in creating Bose-Einstein Condensates (BECs) 
in fermionic clouds. An important step was the develop- 
ment of techniques using magnetically detuned Feshbach 
resonances Jj, which allow to tune the mutual interac- 
tion strength between the fermions over a wide range. 
This novel opportunity to look at a transition from a 
weakly attractive Bardeen-Cooper-Schrieffcr (BCS) state 
to a strongly attractive BEC in one and the same system 
makes it interesting from a many-body point of view (see 
for a recent review). Measurements of the interaction 
strength of a fermionic gas near a Feshbach resonance 
were made by time-of- flight expansion experiments Q. 
The collective excitations showed a strong dependence on 
the coupling strength as was shown experimentally U Q . 
Other experiments observed condensation 0| and the 
spatial correlations @ of the fermionic atom pairs in the 
full crossover regime. Using a spectroscopic technique 
the pairing gap was measured directly [9,]. The remark- 
able result was that the gap values were found in good 
agreement with a simple BCS expression in the whole 
crossover regime. 

One way to access the nature of the many-body state 
is to consider its statistics. A number of works proposed 
to use noise and higher-order correlations to probe the 
many-body states of ultracold atoms [13, [HI, [12, [H, [i3l ■ 
For the BEC-BCS transition the density and spin struc- 
ture factor was calculated Schemes to measure the 
spatial pairing order interferometrically were proposed 
based on correlations in different output channels [l^ . 
To look at pairing fluctuations of trapped Fermi gases has 
been proposed in Experimentally in [18| the spatial 
structure of an atomic cloud was observed directly. This 
enables to determine the density fluctuations for example 
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by repeating the experiment many times or by extracting 
densities at different positions in a homogeneous system 
to obtain the statistics. The shot noise of an atomic beam 
has been experimentally investigated both in bosonic and 
fermionic systems [H, [H [H [li, [H, [H, HI . Further 
aspects of full counting statistics in ultracold atomic sys- 
tems are discussed in the experimental work in Ref. 1261 
and the theoretical papers [27, 28, 29i. [sol. IsH . 

Recently, Amico and coworkers have considered the 
exact solution of the BCS model in some systems us- 
ing the algebraic Bethe- Ansatz [S^] . Explicit expressions 
for average occupations and the number correlators have 
been obtained. Subsequent work has tackled the problem 
numerically and found the Bethe-Ansatz solution to be 
numerically expensive (ssj . The approach to the occupa- 
tion number correlators through the Richardson solution, 
which we develop below, can lead to a numerically less 
expensive method in some cases. A recent review of the 
limit of large particle numbers of the Richardson solution 
can be found in 34. 

In a previous publication, we have calculated the full 
statistics of particle number fluctuations in ultracold 
fermionic gases using a grand-canonical approach [l^ . 
The idea was to consider a 'bin', i.e., a small subsys- 
tem of a homogeneous gas which contains a macroscopic 
number of particles, such that the surrounding atomic 
gas serves as the particle reservoir. Fluctuations can in 
principle be accessed experimentally by performing a se- 
ries of measurements of the number of particles in the 
subsystem at a fixed interaction constant, or by con- 
sidering different bins of the system. The statistics in 
the whole BCS-BEC crossover is hence obtained if one 
sequentially performs such sets of measurements from 
small to large interaction constant. Due to its effec- 
tive single-particle form one can calculate correlation 
functions using the (grand-canonical) Bardeen-Cooper- 
Schricffcr (BCS) ground state solution It was found 
that the BCS-BEC transition yields a crossover in the 
statistics of the particle number. Fluctuations around 
the average particle number are strongly suppressed on 
the BCS side and the statistical distribution is binomial. 
On the BEC side, fluctuations are strongly enhanced and 
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are described by a Poissonian statistics. 

Since real ultracold gases consist of a finite number of 
particles, the grand-canonical approach may be inappro- 
priate. In this article we thus focus on particle-number 
correlations obtained from the exact ground state using 
the methods developed by Richardson [H, 113, IH, HI, 
liol . l4lj | . The Richardson ground state is an eigenstate of 
the total particle number operator N = fii and thus 
allows a canonical treatment of the system. Due to the 
complexity of the Richardson solution, we restrict our in- 
vestigations to a simplified level distribution consisting of 
two multiply degenerate energy levels. This allows us to 
compare explicitly the thermodynamic limit of the exact 
ground state and the approximate BCS solution. Some 
model-independent properties of the statistics can be ob- 
tained in the limiting cases of vanishing or very strong 
interaction. 

Although this is a toy model, it is relevant for a num- 
ber of experimental situations. First, we note that the 
large degree of control possible in ultra-cold atomic sys- 
tems, e.g. by using optical lattices or atomic chips, will 
make it possible to produce few-level systems, which can 
be loaded with a fixed number of fermions in a controlled 
way. Experiments on particle number correlations in such 
systems will be described buy the theory developed be- 
low. Our results apply equally well to level configura- 
tions, in which only two groups of levels are relevant. 
The level spacing within each group has to be smaller 
than the interaction constant; the transition from weak 
to strong coupling appears for interactions of the order 
of the energy spacing between the groups. We would also 
like to mention, that the results obtained below for the 
strongly interacting limit are valid for any level config- 
uration, in which the maximal level spacing is smaller 
than the interaction constant. Hence, we believe that 
systems corresponding to the model we study can be ex- 
perimentally produced, or, at least, some predictions can 
be tested in the limiting case of a strong interaction in 
an arbitrary level configuration. 



tor totally in terms of ^-operators: 
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This is true only in the subspace of fully paired states, to 
which we will restrict ourselves here and in the following. 
The J-th (with J — 1: ground state, J = 2: first excited 
state, etc.) iV-particle eigenstate of ^ has the form 
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Because of (fel)^ = only those terms contribute to the 
sum for which all of the indices /i, . . . , /^v are distinct. 
The coefficient is given by 
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where denotes the sum over all permutations P{i). 
The normalization constant C can be determined apply- 
ing a standard determinant method [37j . The quasi- 

energies -E^"^' in Eq. ([4]) are the solutions of the set of 
coupled root equations 
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In general, the E^^^ are complex quantities, however they 
always appear in complex-conjugate pairs. The corre- 
sponding energy eigenvalue is given by 
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and is thus real, as required. 



The two-level model 



II. PROPERTIES OF THE EXACT SOLUTION 

A. General case 

We start by recalling the basic properties of the 
Richardson solution to the reduced Hamiltonian [H, [13] . 
The Hamiltonian in second quantization and momentum 
space has the form 
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where &/ = c/|C/| are hard-core bosonic annihilation op- 
erators. The reduced Hamiltonian captures only the scat- 
tering fermions which occur in time-reversed states. It is 
therefore possible to express the particle number opera- 



We will now consider the special configuration involv- 
ing A'^ particles in two multiply-degenerate energy lev- 
els eo and t\ with the degeneracies fio and fii, respec- 
tively . In the following the subscripts and 1 will re- 
fer to one of the lower levels and one of the upper lev- 
els, respectively. The coefficients Lp from Eq. ^ can 
be expressed in terms of a single variable v that indi- 
cates the amount of particles in the upper level ei [s^ : 
. . . , /jv) — > (/^(•^^(i^). This leads to a simplifica- 
tion in finding the Richardson solution. Introducing the 
abbreviations 

^ 2iVeo + 2i/(ei - eo) - .gi^(f^i - I/+ 1) 

-g(iV-z/)(r!o-iV + /^+l) , (7) 

A, = g{N -v){^^-v) , (8) 

= gv{^Q - iV + I/) , (9) 
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the coefncients are determined by the continued-fraction 
formula 
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(10) 

ip^'-'' (0) has to be extracted from the normahzation con- 
dition 

i:(/-.)C;)^'"<"'i'^i- (") 



j/=0 



One can find an expression for the total energy ^ ap- 
pearing in Eq. pop directly from the root equation 
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without having to resort to the quasi-energies from 
Eqs. ©. 

In the following, we will use the level spacing ei — eo = 1 
as our energy scale. The only (dimensionless) parameter 
left to characterize the system is the ratio between the 
interaction constant G and the level spacing. 

The model of two highly-degenerate levels is clearly 
an artificial model, which cannot be mapped to generic 
many-particle systems. However, we believe the model 
is nonetheless relevant also for experimental realizations 
due to two reasons. First, we show below that devia- 
tions of the occupation number correlators from a sim- 
ple BCS mean-field treatment are relevant and can be 
detected in not too large systems. In experimental real- 
izations of interacting Fermi systems in ultracold atomic 
gases an unprecedented variability of system parameters 
has been experimentally demonstrated 2]. We therefore 
hope that our investigation will stimulate experimental 
efforts to create few-particle strongly-interacting systems 
and study the effects we predict below. Second, within 
the same reasoning we believe that the large variability 
of tailoring atom systems in (magneto-)optical traps or 
via atomic chips will make it possible to create an arti- 
ficial highly-degenerate two-level system and to use it to 
study in a controlled manner the transition to the ther- 
modynamic limit in a particularly simple system, as we 
predict here. 



C. Example 

We want to illustrate some characteristics of the 
Richardson solution by means of a simple setup within 
the two-level model. We consider, therefore, a system of 
two hard-core bosons in two threefold degenerate levels 
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FIG. 1: (Color online) Development of the ground state and 
first and second excited state energy vs. interaction constant 
g for a simple two-boson system. The two levels at eg = 
and ei = 1 are both threefold degenerate. 



(r^o = r^i = 3) of energy eo = and ei = 1. Figure [T] 
shows the three root solutions of Eq. as a function of 
the interaction g. Obviously, in the non- interacting limit 
at g ^ 1, the energies reduce to the bare pair energies 0, 
2 and 4, corresponding to the case that v = Q, v ^ 1 and 
V — 2 particles respectively are in the upper level. With 
increasing interaction g « 1, the ground state and first 

excited state energies are lowered continuously, whereas 

(3) 

the second excitation energy approaches £3 2 and is 
then independent of the interaction constant g. 

In the following, we concentrate our investigations on 
the behavior of the ground state. Figure O shows the be- 
havior of the many-body occupation [(^s^^-* (i^)] as a func- 
tion of g. At vanishing interaction only the lower three 
energy levels at eo are occupied. From the normaliza- 
tion condition ^ it thus follows that [ip^^\Q)f = 1/3, 
since there are (2) = 3 distinct possibilities to distribute 
two particles among three levels. The average occupa- 
tion number of a lower level is hence given by (rio) = f ■ 
At strong interactions g 3> 1, all levels tend to become 
equally occupied. In this limit, we can therefore neglect 
the level spacing and consider simply a single energy level 
with a total a single total degeneracy Q. = Vt^ + — Q. 
Equation ^ simplifies to [(^(1) (i/)]~2 ^ (-n^ ^ Q ^ 



which is independent of v, and the average particle num- 
ber of a level is given by (no) = (ni) = N/Vl = 1/3. We 
want to point out here that the equal occupation of all 
levels in the strong interacting limit is not restricted to 
this specific level model but rather a general feature of 
the Richardson ground state. 

For g > 1 it follows from Eq. ^ and Eq. ^ that 
the ground state energy approaches 
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2iVeo - gN{n - N + 1) , (13) 



V^^Ho) 

since if (i)(l) « (fi^^HO). This useful relation can e. g. be 
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FIG. 2: Many-body coefficient [ip'^^\u)]'^ as a function of n 
for the ground state energy e'^\ 

taken as an initial energy guess in the whole interaction 
regime, when it conies to finding the roots of Eq. (fT^ . 

III. PARTICLE-NUMBER CORRELATIONS: 
GRAND-CANONICAL VS. CANONICAL 

We address now the correlations following from the 
exact ground state. We particularly focus on the differ- 
ences that occur in a canonical treatment compared to 
applying the grand-canonical BCS solution. At first, we 
define the particle number cross-correlator between the 
occupations of levels / ^ /' 

gifJ') ((^/-("/))(^/' -("/'))) 

= {nffif) - {nf){nf,) , (14) 

which represents a direct measure of how much the par- 
ticle number of level / fluctuates around its mean value 
in the presence of a fluctuation around the mean value 
of the particle number of level /'. 

The grand-canonical BCS wavefunction [35] is given by 

\BCS) =l[{uf +vfb\)\0) , (15) 
/ 

with vj^l- u) = (1 - {ej - M)/v/(e/-M)2 + A2)/2, 
where ^ is the chemical potential. The mean field A and 
^ are fixed by the self-consistency equations 

l^^-gY^UfVf , N^2Y,v} . (16) 
/ / 

The simplification by the mean field Ansatz, that is the 
reduction of the many-body interaction to an effective 
one-body interaction, has a direct consequence on cross- 
correlations: Since (fT5|) is a product state the different 
level occupations are uncorrelated and, hence, g{ f^ /') = 
0. As we will see in the following sections, the correlations 



will be non-zero if the many-body interaction is taken 
into account beyond the mean-field approach. 

Due to the operator identity fif = for Z = 1,2, ... 
in the subspace of paired particles, the auto-correlation 
function of a level, Eq. (|14p with f — f, is totally de- 
termined by its average particle number and thus does 
not contain any additional information. In the following, 
we will therefore concentrate on the investigation of ex- 
act average particle numbers and exact particle number 
cross-correlators in the form of Eq. (|14p . 

A. Exact correlators in the two-level model 

We now determine the explicit form of the particle 
number cross-correlator (|14p in the two-level model. We 
only have to consider three different kinds of correlators, 
since all degenerate levels are equivalent. If two levels of 
the same energy are distinct, we indicate this by priming 
one of the indices labeling the energy of the level. The 
three different cases take the form 

5(0,0') = {nono')-{nof (17) 
3(0,1) = (noni) - {na){ni) (18) 
5(1,1') = {h,h,,)-{h,f (19) 

with (assuming that N < Hq, Hi) 

<^».-=.>^i(/_\\^)r;:;)h'"<"f. 

and 

If we allow particle numbers exceeding one or both 
degeneracies, the boundaries in the sums appearing in 
Eqs. pT|) . (PU)) and (PT|) have to be adjusted accordingly. 

B. Relation to counting statistics 

In the two-level model, we are also able to spec- 
ify the full statistics of the occupation numbers. This 
quantity can in principle be obtained by measuring re- 
peatedly the occupation numbers and finding the prob- 
ability P{nf,nf') that two levels / and /' have oc- 
cupations n f and n f . This full statistics can be 
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equivalently expressed through the cumulant generat- 
ing function Sff,{xf,Xf') = lnE„,,n^., exp(ix/"/ + 

*X/'"/')-P//' ("■/)"/') Hi- The cumulant generating 
function for hard-core bosons in a fully paired state is 
given as a function of two counting fields x / and x f '■ 

gSff,{xf,Xf') ^ ^gi(x/ft/+x/'ft/')) (23) 

= 1 + (hf) {e^^f - 1) 

+ {nf') {e'^'f' - 1) 

+ {nfnr) (e^« " l) (e^^^' - l) ■ 

Consequently, the only correlator which needs to be 
known to fully determine the CGF is the one in the last 
line of Eq. ((23)) . for which we are able to give explicit 
expressions here, due to the simplicity of the model. 

In the case of non-interacting particles, e.g. hard-core 
bosons in the BCS mean- field treatment, Eq. factor- 
izes according to 

f,Sff,{xf,Xf') ^ gS/(x/)+Sj./(x/') (24) 

= [1 + (ny) (e'« - 1)] 

X [1 + (uf,) {e^^f - 1)] . 

This is the CGF of uncorrelated particle numbers. Com- 
paring these general results for the counting statistics 
with the correlators discussed in the previous subsection 
we observe that in this special case the counting statis- 
tics contains not more information than the correlators 
alone. Or, in other words, if the correlators Eq. (fT7)) - 
are known, one can use Eq. to calculate the full 
counting statistics. 
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FIG. 3: (Color online) Average particle number of one of the 
upper levels as a function of G for various system sizes. The 
black solid lines correspond to the solutions obtained from 
BCS theory. Note that the self-consistency Eqs. (|16|l in this 
case only have real solutions for G > 1. The average particle 
number of one of the lower levels follows from (no) = 1 — {hi ) . 



that all coefficients [tp^^^ (/i • ■ ■ Jn)] become equal in the 
strongly- interacting limit, (j26p is a universal property of 
particle- number correlators, which is valid also for arbi- 
trary level configurations and not only restricted to this 
simple model. 



D. The two- level model at half filling 



C. Asymptotic behavior 

Before we discuss the general results for an arbitrary 
interaction constant, we obtain analytical expressions for 
the correlators in the limiting cases of weak and strong 
interactions. This is possible since the coefficients (|4]) can 
be directly determined from the normalization condition 
(fTTjl without having to solve the root equation, Eq. (H^). 
In the following, we will assume that N < flo, ili. For 

g < I, Eq. (ini) reduces to [^(i)(0)]"^ = (^^°). From 
Eq. (HOD and Eq. ^ thus follows 



ff(0,0') 



N - (no) ' 



(25) 



with the system-size-independent average particle num- 
ber (no) — N/Qq. The remaining correlators are zero. 

For g ;» 1 we have correspondingly [(^(^^(i/)] = (^), 
where 51 = f^o -I- ^^i. Hence we get 



in) 



N - (n) 



(26) 



Again (h) = N/D, is the system-size and interaction- 
constant independent average occupation number of ev- 
ery level. Since it is a feature of the Richardson solution. 



We now discuss the numerical results of the aver- 
age particle numbers and correlators above as we ap- 
proach the thermodynamic limit starting from finite sys- 
tem sizes. In the evaluation of e. g. Eq. (PO)) and Eq. (|2ip . 
we hence have to assure that the involved quantities scale 
in the correct manner. The continuum limit is obtained 
by taking 17 oo, while leaving and G — gfl constant 
[SI, S, . We call G the 'system-size-independent cou- 
pling constant'. Under these assumptions, increasing the 
particle number will lead to the BCS results in the ther- 
modynamic limit. 

At first, we investigate the two-level model at half 
filling with equal degeneracies of both energy levels, 
viz. r^o = = TV. Figure [3] shows the average parti- 
cle number in one of the upper levels as a function of G 
for various system sizes. We can see that there is already 
a fairly good agreement to the BCS results in the case of 
only 32 particles. Due to the particle-hole symmetry of 
the system, the connection between the average particle 
number of a lower level and an upper level is given by 
(rio) = 1 — (n-i). The average particle numbers in the 
limits of weak and very strong interactions are system- 
size independent: for G ^ 1, only the lower energy band 
is occupied. For G ^ 1 as mentioned in Sec. Ill CI we 
obtain an equal occupation of all levels. 
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FIG. 4: (Color online) ff(0,0') = 1') (upper plot) and 
g{0, 1) (lower plot) as a function of G for various system sizes. 
The inset shows that, in contrast to 17(0,0'), g{0, 1) is always 
negative in the low-interaction regime. 



In Fig. 3] the corresponding correlations are given as a 
function of G. Note that, due to particle-hole symmetry 
in the half-filled case 17(0, 0') = g{l, 1'). The behavior of 
the average particle numbers in the strongly-interacting 
case has a direct influence on the correlations causing 
(7(0,0'), (7(0,1) and (7(1,1') to become equal in magni- 
tude for a fixed system size. At vanishing interaction, the 
lacking possibility of reshuffling particles in a fully occu- 
pied band leads to zero-correlation. A comparison of the 
plots in Fig. 3] shows that for a given TV, the crossover 
happens over a smaller range of G in the case of g{0, 1). 
Obviously, the fact that all occupations 1/ contribute to 
this correlator - contrary to Eq. (fT7|) , where v = N only 
enters through the normalization pT|) - leads to a faster 
saturation with increasing coupling constant. 

As a general feature, one finds that the particle number 
correlators of distinct levels tend to converge to the zero- 
correlation line of the mean-field approach in the whole 
interaction regime as one increases the number of par- 
ticles. A direct indication of the fermionic origin of the 
hard-core bosons is that, at first sight, in the non-limiting 
cases N 00 and G 7^ all correlators are negative, cor- 
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FIG. 5: (Color online) Zoom into the region of positive cor- 
relations of (7(0,0') — g{l, 1') in the weak-coupling range for 
different system sizes. A positive peak develops around G — 1 
for small systems and becomes maximal for « 16. For large 
system sizes the absolute values become smaller again, but the 
overall feature sharpens. 



responding to anti-correlated particle numbers: Due to 
the presence of a particle in level /, it is less probable to 
find another particle at the same time in level /' than in 
the uncorrelated case. However for (7(0,0'), we observe 
a range at intermediate interactions, where particles of 
a certain energy promote other particles to occupy the 
same level. It also leads to another point of vanishing 
correlation for G 7^ and fixed N; see Fig. [Sj Evalu- 
ating (7(0, 0') in second-order perturbation theory shows 
that this effect starts to occur for iV > 2. There is a 
resonance effect with a maximum peak value in the pos- 
itive correlation between 16 and 32 particles. We do not 
find a positive range for g{0, 1). This surprising finding 
is confirmed analytically by a perturbative calculation in 
the appendix. 



E. The two-level model away from half filling 

It is also interesting to study the system away from half 
filling. As an example we now look at the case of quarter 
filling. Again, we assume that Qq = Qi and to have a 
direct comparison to the model at half filling, we chose 
the same system sizes as in the last example. Because of 
particle-hole symmetry we have the following relations 
between filling factors 6 — N/{ilo + fli) and 1 — 6. For 
the average particle numbers 



{ho)s = 1 - (711)1-5 
("-1)5 = 1 - ("0)1-5 



and for the correlator 



ff(0,0')5 = 5(l,l')l-5 

g{l,V)s - .g(0,0')i_5 
5(0,1)5 = 5(0,1)1-5 . 



(27) 
(28) 

(29) 
(30) 
(31) 
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FIG. 6: (Color online) Average particle number of one of the 
upper levels as a function of G for various system sizes in 
the quarter-filled case 5 = 1/4. The BCS Eqs. ([T6|l have real 
solutions for G > in this case. The average particle number 
of one of the lower levels follows from (no) = 0.5 — (ni). 



In the following, we will consider the case oi S = 1/4 
(which is therefore equivalent to 6 = 3/4). The average 
occupation of one of the upper levels for this case is shown 
in Fig. [51 The occupation of one of the lower levels is not 
shown, since it follows from (fio) = 25— {fii). We see that 
the saturation at large interaction constant happens at 
larger G than in the half-filled case (c. /. Fig. [3]). Note 
that here the BCS solution always exists and is indistin- 
guishable from the exact solution already for 16 bosons. 

Figure [7] shows that all correlators are now negative 
and different from each other and approach their lim- 
its in the strongly interacting case more slowly than in 
the half-filled case. f;(0, 0') shows an interesting behavior 
for vanishing interaction that is caused by the partially 
occupied lower energy band allowing particles to change 
states among the lower levels. This leads to a finite value 
also for G ~ and a decay of the correlator with increas- 
ing system size, see Eq. (PS)) . It is also remarkable that 
g{Q,0') is suppressed by increasing the interaction. Also, 
in agreement with the perturbative results Eq. (|A.1|1 the 
effect of the interaction is of second order in the inter- 
action constant. The other correlators show a similar 
behavior as in the half-filled case. 
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FIG. 7: (Color online) g{0, 0'), ff(0, 1), and 3(1, 1') as a func- 
tion of G for various system sizes at quarter-filling 5 — 1/4. 



IV. CONCLUSION 

We have investigated exact particle-number correla- 
tions of ultracold fermionic gases in a canonical Ansatz 
using the Richardson solution. By means of a special 
configuration involving two degenerate energy levels, cor- 
relation functions have been derived and evaluated nu- 
merically for different mutual interactions between the 
atoms and different system sizes. The particle numbers 
in different levels turn out to be mostly anti-correlated 
revealing the fermionic origin of the paired particles (the 



hard-core boson property). Approaching the thermody- 
namic limit, those correlators decay to zero in the whole 
interaction regime. This is in agreement with the predic- 
tions of BCS theory. In the limit of strong interactions we 
were able to give closed expressions for the correlations, 
which are also valid for the general case of arbitrary level 
configurations. Due to the complex algebraic structure 
of the Richardson solution, only a comparatively special 
model could be investigated in this work. The discus- 
sion of more general systems remains an open problem. 
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Nevertheless, we believe our predictions can be tested in 
tailored few-particle systems of interacting fermions, e. g. 
with atomic chips. 
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The Richardson solution and correlators can be found 
perturbatively in the interaction constant. For N < 
r^o, ^1 we find the expression 



5(0,0') 



iV(»o - N) 



_ GV l]i(»o-jV+l) 
N 2 



Nni{no-N + i) /G 
i6rjg(i7o - 1) 



2N{N - 1){5N - 6) - {7N - 2f7o - 7)(3iV - - 2)no 



+ {13N{N - 1) + 7f7o(2 - 3iV + r2o))r2i + 2{N - f7o)f^? 



(A.l) 



For a half- filled band {N = = fli) the zeroth and the second-order term vanish and the expansion to 6th order 
yields 



5(0,0') «( 



1 



1 



167V2 



8iV3 



13 



167V2 327V3 i6iV4 



+ (D{G^) 



(A.2) 



The G'* term is positive for N > 2, but gets smaller for increasing N. For increasing G the 6th-order term takes over 
and leads to a negative correlator in the end. 
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